Enhanced shear separation for chiral magnetic colloidal aggregates 
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We study the designing principles of the simplest colloidal propeller, an architecture built from 
four identical spheres that can couple translation with rotation to produce controlled drift motion. 
By considering superparamagnetic beads, we show that the simultaneous action of a magnetic field 
and a shear flow leads to the migration of the cluster in the vorticity direction. We investigate the 
dependence of the migration velocity on the geometrical parameters of the cluster, and find that 
significant cluster separation can be achieved under the typical operation conditions of microfluidic 
devices. 

PACS numbers: 82.70.Dd; 47.57.cf; 85.70.Rp 

The motion and the propulsion of microscopic living 
organisms, the development of top-down techniques to 
carve micrometer-size devices like pumps, switches and 
motors, the supramolccular chemistry quest for molec- 
ular motors and the new bottom-up techniques that 
emerge from the Soft Nanosciences, have lead to a reborn 
interest for the hydrodynamics at low Reynolds numbers 
Re. In the microscopic world of fluids where typical di- 
mensions £ are of the order of, or smaller than a few mi- 
crometers, I % 10 /im, where the velocities v are smaller 
than a few micrometers per second, v < 10 /im.s" 1 , the 
ratio Re = £vp/r) between inertial and viscous forces is 
extremely small, Re < 10 -4 . At these scales, inertia is 
irrelevant: for all practical purposes when the applied 
forces stop, the movement stops. Designing feasible and 
efficient propulsion devices for low Reynolds numbers is 
arguably one of the main current scientific challenges 
in the Nanosciences. Inspired by recent progress in the 
self- assembling small numbers of colloidal particles into 
aggregates with well defined geometries [l|, and moti- 
vated by the possibility of manipulating colloidal aggre- 
gates of superparamagnetic beads with external magnetic 
fields 0], we propose and study in this Communication 
the simplest feasible colloidal aggregate that will behave 
as a small propeller, capable of transducing rotation into 
translation and vice-versa. 

The propeller effect consists in rotating at constant 
angular velocity uj a chiral aggregate around a fixed ori- 
entation axis, conventionally named z, and subsequently 
referred as the vertical axis, causing a net constant drift 
velocity v z along z. The propeller efficiency of an aggre- 
gate of size L can be expressed as a dimensionless pitch 
P p = 2itv z /(Luj) that, in the Stokes flow approximation, 
provides a size independent characterization of its ability 
to couple rotational and translational motions. Spatial 
preorientation of the aggregate is required to observe a 
sizeable linear propeller effect, otherwise strongly sup- 
pressed by the free rotational Brownian motion. We show 
below that this can be achieved with a constant uniform 



magnetic field acting on superparamagnetic colloids. 

Rotation of the aggregate can, in principle, be enforced 
by means of circular polarized light Alternatively, a 
precessing magnetic field with small but finite tilt angle 
/3 could also be used, but only with low efficiency (i.e.^ 
/3 2 ) In any case, it is simpler to rely on flows to rotate 
the clusters, leading naturally to a shear induced vertical 
drift in Couette geometry. 

It was recognized some time ago that the migration 
under shear flow of a colloidal cluster with randomized 
orientation only shows a cubic dependence in the shear 
rate 7 [B| . Shear separation becomes linear in 7 as soon as 
aggregates get prcoricntcd with respect to the shear flow, 
with z parallel to the vorticity direction. We introduce 
the shear separation efficiency ratio P s = —At:v z /{L^i) as 
a counterpart of P p , based on a correspondence u = —7/2 
between angular velocity and shear rate. Propeller and 
shear separation are related, but were found incquivalent 
in the cases investigated so far, and comparing P p and 
P s is an important issue. Strong constant shear rates can 
be experimentally applied, e.g. using soft lithography 
microfluidic devices [g, Q and shear separation has the 
potential to become a method of choice for partitioning 
aggregates according to their chirality. 

We first discuss the preorientation of the colloid aggre- 
gates by the applied field. Each superparamagnetic bead 
i making up the aggregate acquires an induced magnetic 
moment rhi and the equilibrium position of the aggregate 
is the one minimizing the sum — J^. m, • B/2, where B 
stands for the uniform, gradientless, magnetic induction 
in the absence of colloids. We also assume that the mag- 
netic content of the beads is homogeneous, and that the 
aggregate is rigid. The actual determination of the best 
orientation of the aggregate depends on the strength of 
mutual magnetic moments interaction, controlled by a 
dimensionless parameter e = /ioa;>/(327ra 3 ), with po the 
permeability, on, the magnetic polarizability, i. e the ratio 
between magnetization m and induction B, and a the ra- 
dius of the beads. Provided that i) the applied magnetic 
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field ||-B|| is less than a critical value B c and does not 
bring m, close to its saturation value, ii) the time varia- 
tions of B are slower than the Neel relaxation time r^r of 
the magnetic particles and Hi) the interaction parameter 
s remains smaller than unity, then the magnetic energy 
reduces to E m (B, ft) = J27<j wu " n a dipole-dipole in- 
teraction equal to Vij = 8aba 3 eB 2 r~j 3 (l — 3 cos 2 % ) , 9ij 
the angle between the magnetic field and the vector re- 
connecting beads i and j, and Q a set of Eulcr angles 
fixing the orientation of the cluster. This dipolar en- 
ergy is the leading term of an expansion in powers of e 
that fully takes into account the magnetic interactions 
between beads. 

The quadratic dependence of E m (B,Q) in B at fixed 
f2 leads to a unique optimal orientation of the aggre- 
gate with respect to B if the spectrum of the quadratic 
form has a single lowest eigenvalue, which is the case 
for generic, non-symmetric aggregates. Otherwise there 
may be a degenerate subspace of low energy configura- 
tions, while aggregates with high symmetry end up being 
totally insensitive to the applied field, as it happens e.g. 
for a regular four beads tetrahedron. 

For dimers, the preferred orientation is along the field 
direction; for trimers, as for any planar aggregate, the 
plane of symmetry of the cluster contains the field di- 
rection. These are particular instances of more general 
orientation rules which are strongly constrained by the 
symmetry properties of the aggregate. Assuming that a 
given aggregate possesses a well defined orientation, then 
any symmetry plane a of the aggregate is either normal 
to B or contains it, any rotation axis C n of order n > 3 
is parallel to B, and any rotation axis Ci ends up being 
parallel or orthogonal to B. We assume from now on 
that our colloidal aggregates possess a well defined ori- 
entation axis and rotate around it. In the case of shear 
induced migration, the vorticity of the shear flow should 
bc aligned parallel with B in order to facilitate the rota- 
tion of the aggregate. Flows with vorticity perpendicular 
to B induce torques working against the magnetic stabil- 
ity of the aggregate, ending up with no or little rotation, 
and are not considered further in the present approach. 

Let us consider the linear relation v z = kuj between 
drift and angular velocities of an aggregate with con- 
strained rotation axis z. If an inversion center of the 
aggregate exists, this inversion symmetry preserve both 
k and uj while reversing v z . If a mirror symmetry plane 
containing the axis z exists, it preserves k and v z . but 
reverses uj. If the mirror plane is orthogonal to the z 
direction, it preserves k and uj, and reverses v z . There- 
fore, k vanishes for all non-chiral aggregates. The same 
consideration holds when a torque T is applied instead 
of a constant angular velocity uj. Dimers and trimers 
are de facto ruled out as they possess a mirror plane 
containing the field. It requires at least four beads to 
build a propelling aggregate, the only point-symmetry 
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FIG. 1: Sketch of the colloidal propeller, made of identical 
superparamagnetic beads. For displaying purposes only, the 
beads are represented by spheres of a diameter smaller that 
the arm length in (b) and (c). The propeller consists of a 
main axis and two perpendicular arms that form an angle 8 
(a). For small values of S, one stable equilibrium position 
under vertical magnetic field, that we name A orientation, is 
shown in (b), the other being obtained by a rotation of 180° 
around the main axis. For larger values of 5 the equilibrium 
position that we name £ orientation is shown in (c). 



groups allowed for propelling aggregates are, according 
to Schoenflies terminology, Ci, C„ and D„ (n > 2) Q. 

We will now show that some tetrameric assemblies dis- 
play indeed propeller and shear separation properties, 
and compute the magnitude of these effects. We consider 
an aggregate of four identical beads rigidly connected, as 
depicted in Fig. ((TJa). In this configuration, bead one 
and bead two at a diameter distance L = 2a define the 
main axis. The third bead, connected to bead one defines 
a lower arm of length L perpendicular to the main axis, 
and the fourth bead, connected to bead number two de- 
fines an upper arm of length L, also perpendicular to the 
main axis. The upper and lower arms are separated by 
an angle S. Aggregates with angles S and —S only differ 
in their chirality, also called enantiomers. We consider L 
as the characteristic size of the aggregate. 

The equilibrium position of such cluster under mag- 
netic field depends on 5. For small values 5 < S c , the 
typical orientation of the cluster that we name orienta- 
tion A is displayed in Fig. (jTJb) . One observes that the 
field direction is along the (unique) axis of symmetry Ci 
of the tetramer. For larger values of 6, the stable con- 
figuration, that we call orientation E, adopts an helical 
conformation as shown in Fig. ((He), with B perpendic- 
ular to the C2 axis, in accordance with the above ori- 
entation rules. The A-E stability exchange, caused by 
an eigenvalues crossing of E m (B, fi) takes place around 
S c ~ 66.4°. 

Within the Rotne-Prager-Yamakawa dipolar approxi- 
mation Q the velocity of the particle located at R m is 
given by the hydrodynamic relation V r „ = J2 n H TO « 'F n , 
where F„ are the forces acting on particles and the mo- 
bility matrix H m „ is approximated by H„„ = 1/(671770) 
and H„ m = H(R m — R„) for m ^ n with H (r) = 
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l/(87T77r 3 ) (r 2 l + rr + 2a 2 (l/3-rr/r 2 )). The validity 
of such modeling of cluster hydrodynamics was investi- 
gated in and suggests an accuracy within 30% for 
the rotational drag coefficient. The rigid cluster move- 
ment is reducible to three translational and three ro- 
tational degrees of freedom collectively denoted as Q a , 
a = x,y,z,4>,6,il>, to which are associated generalized 
forces and torques T a . Internal constraints can be elim- 
inated by suitably projecting the inverse hydrodynamic 
tensor H _1 over the collective displacements dTL m /dQ a , 
leading to 0: 



Qa — 



7=1, ...6 



h (6) T 



(1) 



where hi 6 7 is a 6 x 6 symmetric mobility matrix re- 
stricted to the collective coordinates obeying (h _1 ) afc = 
J2 m .n dRn/dQa • (H _1 ) nm • dR m /dQ b , F 1 arc the exter- 
nal forces and torques components, and lA a the general- 
ized components of the streaming flow. In the ideal pro- 
peller case, the motion reduces to only two coordinates 
9 and z, standing respectively for the rotation around z 
axis, and translation parallel to it. In the absence of flow, 
one gets: 



K a 16, 



h^F, 



(2) 



with F z the force, Te the torque and h^ 2 ) a 2 x 2 sym- 
metric matrix. Eq. ([2]) suffices to discuss the properties 
of propeller configurations, while cq. ([T]) must be numeri- 
cally integrated when considering shear separation exper- 
iments or rotation under precessing magnetic field. Nu- 
merical integration of these equations by a fourth-order 
Runge-Kutta method [l2| provides complete solutions for 
the time evolution of the center of mass and orientation 
of the propeller. 

Applying a torque 7g and no force F z results in an he- 
lical trajectory 0(t), z(t) associated with a pitch P p = 

(2) (2) 

(2n/L)z/0 = (2Tr/L)hy'Jh)f> g which confir ms the geo- 
metrical, time-reparametrization invariant of the motion. 
In a Couette flow with vorticity parallel to the magnetic 
field, the aggregate rotates producing an oscillatory verti- 
cal movement with a net macroscopic migration along the 
vorticity direction. The angular velocity of the propeller 
also oscillates in time, as shown in Fig. ©, depending 
on the relative attitude of the cluster with respect to the 
anisotropic shear flow. Curiously enough, the faster ver- 
tical motion coincide with the slowest angular velocity 
part of the trajectory. This can be understood if one ob- 
serve that an anisotropic shear flow can exert a lift force 
on a non rotating aggregate. 

We show P p as a function of S in Fig. (|3|). For the 
two points with no chirality (6 = 0° and 180°) the pitch 
vanishes as expected. P p is also seen numerically to van- 
ish around 6 a ~ 95°, thus separating the E regime in 
two regions with pitches of opposite sign. A maximum 



0.4 



0.3 



0.2 



0.1 



0.0 



-Z Q) 
&_ 

0) 

-3 S 



-4 S 

CD 
T3 



-90 -180 -270 
9 (degrees) 



-360 



FIG. 2: (Color online) Angular velocity (red dashed line) and 
vertical position of the center of mass (black full line) of the 
propeller as a function of time while it turns in a shear flow. 
The straight dotted line represent the vertical motion of a 
propeller rotating with same angular velocity, S = 40°. 
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FIG. 3: (Color online) Pitches P s (shear flow, full black line) 
and P p (propeller, dashed red curve) as a function of the 
angle S for touching beads L — 2a. A and E orientations 
are depicted on the graph, separated by a critical value 8 C . 
Inset: stability A 9 of the magnetically oriented aggregate. 



of P„ 



p — 0.05 corresponds to S c 
of P v ~ —0.05 obtained for S m 



66° and a minimum 
130°. As a complete 
rotation of the cluster corresponds to a vertical shift of 
Az = P p L, it requires about 20 rotations of our cluster 
to move upwards over its own size. 

The shear flow pitch P s correlates well with P p in the 
A orientation, but has a sign opposite to P p in the E ori- 
entation. In particular, P s and P p apparently vanish for 
the same value 5d, in agreement with a general argument 
stating that P p = => P s = 0, without reciprocity Q. 
The propeller efficiency in the absence of external flow 
is overall lesser than the shear slow efficiency in A orien- 
tation while it is stronger in E orientation. The relative 
sign of the two effects depends on the details of the col- 
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loidal cluster. 

The inset of Fig. |3] shows the stability X g (d) of the 
magnetic field induced orientation of the cluster, as a 
function of 5. X g is defined as twice the difference be- 
tween the two smallest eigenvalues of the quadratic form 
E m (B,Q), with 8a e = 1. As a result, the depth of 
the magnetic induced energy minimum and the maximal 
equilibrium restoring magnetic torque can both be ex- 
pressed as E mag = Tmag = ^ea b X g B 2 . 

We have so far neglected thermal forces in our anal- 
ysis. The competition between drift and diffusion de- 
fines a Peclet number Pe = \v z \L/D t , with a diffusion 
coefficient D t = ksT / (QirqL). In the sheared situation 
v z = ~P S ^L 1 47r and the ratio between vertical and Brow- 
nian displacements scale as (jPePst/Sn) 1 / 2 , defining a 
characteristic time t c = 87r/(7PeP s ) beyond which the 
drift displacement dominates. 

Let us illustrate with numbers these phenomena in 
the case of a shear experiment. Many common su- 
perparamagnetic beads contain crystallites of 7FC2O3 
with estimated radius r c ~ 10 nm (l3j . Each one 
of these magnetic grains bears a thermally fluctuating 



moment of magnitude m p 



1.9 x 10" 18 A.: 



with 



random orientation. This sets the saturation field to 
B c = ksT /m p ~ 2.2 mT and the linear response polariz- 
ability of crystallites amount to a p ~ 3 x 10 -16 A.m 2 .T _1 . 
The Nccl relaxation time was experimentally found close 
to tn ~ 4 x 10~ 5 s The polarizability a b and cou- 

pling parameter e of a bead both depend on the mag- 
netic material volume fraction leading respectively to 
a b = a p (a 3 /r 3 )$ and e ~ 3.6$. 

The dipolar energy used in our approach is correct for 
the moderately magnetic beads (<& = 0.02, e ~ 0.07) 
used in 15j while corrections to the magnetic energy are 
required for the stronger magnetic beads ($ ~ 0.6, £ ~ 
2.2) studied in [14j. The magnetic orientation energy 
must excede the thermal energy ksT. Our calcula- 
tion is restricted to the linear polarization regime, lead- 
ing to a magnetic stabilizing energy of order £ mag ~ 
4.8 x 10 6 a 3 $ 2 A s fcsr, with a expressed in micrometers. 
The magnetic torque Tmag must match the viscous torque 
Tdra g ~ 4irr]L 3A f exerted by the flow, with a dynamic 
viscosity r\ = 10 _3 Pa.s in the case of dilute water solu- 
tions. This, at room temperature, limits the shear rate 
to ^max ~ 2.0 x 10 5 $ 2 X g irrespective of cluster size L. 
The larger drift velocity v z consistent with these con- 
straints scales as V max — P s L ^max/ (47r), or equivalently 
V max — 3.1 x 10 4 P s a$ 2 A g , with a corresponding Pe scal- 
ing as Pc = 5.6 x 10 5 P5a 3 $ 2 A 9 (a in micrometers). 

For micron sized beads with moderate content in mag- 
netic material ($ = 0.1), assembled with an angle 6 = 
40°, one reads from Fig. [3] the values A 9 ~ 1.4 and P s = 
0.25. One obtains a maximal shear rate j m ax ~ 2800 s _1 , 
a maximal drift velocity V max = 110 fim.s^ 1 , t c = 18 /xs 
and Pc=2000. This should be fast enough to induce some 
cnantiomer separation using a microfluidic chip like the 



one described in [6j with a magnetic field aligned parallel 
to the chip and orthogonal to the main direction of the 
flow. The scaling Pe ~ a~ 3 limits our results to sizes 
a = 80 nm below which Pe < 1 and Brownian forces 
cannot be neglected any longer. 

In summary, we have applied an hydrodynamic inter- 
action formalism for studying the rotational-translational 
coupling for a chiral propeller made of four identical 
beads. This coupling can be optimized by changing the 
geometry of the aggregate. Chiral superparamagnetic 
colloidal enantiomers can be oriented and separated un- 
der the simultaneous action of a magnetic field and a 
shear flow even at relatively low shear rates since they mi- 
grate in opposite directions with velocity proportional to 
the applied shear. Propeller and shear separation prop- 
erties are distinct features, their relative sign may change 
depending on the details of the cluster, and are not sim- 
ply related to the handedness of the chiral aggregate [l(| . 
Quantitative estimates of the magnitude of this magnetic 
field assisted drift are given that support the feasibility 
of the experiment with microfluidic devices. 
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